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Abstract 

fS We extract the strong coupling constant a s (m 2 ) from the semileptonic r decay 



Oh 



width taking into account resummation effects from the running of the strong cou- 



pling constant. In the MS scheme the result reads a s = 0.375 ±0.007 to third order 
and a s = 0.378 ±0.007 to fourth order in the /3-function, respectively, where we use 
the recently computed four- loop coefficient fa. These values for the coupling con- 
stant have to be compared with the value a s = 0.354 ± 0.005 derived from a third 
order analysis of r decays. We determine the exact value of the convergence ra- 
dius of the perturbation series by analyzing the singularity structure of the complex 
coupling constant plane. 



The precise determination of the strength of the strong interactions is of paramount 
importance for phenomenological applications of the Standard Model. Within perturba- 
tive QCD the strength of the strong interaction is quantitatively measured by the coupling 
constant a s . Although the value of a s depends on the choice of renormalization scheme 
in finite orders of perturbation theory, the running coupling constant a s (Q 2 ) has now 
become a widely recognized measure of the strength of the strong interactions QTJ] . 

It has been realized since many years that r decays provide an ideal setting for obtain- 
ing accurate information on strong interactions in the low energy domain. The inclusive 
semileptonic r decay width was analyzed in detail within perturbative QCD including 
also the small nonperturbative power suppressed effects. A numerical value of a s was 
obtained from this analysis [[|. The error that was quoted in these analyses, however, 
accounts mainly for the statistical uncertainty in the experimental decay width T T and 
does not include systematic errors due to the truncation of the perturbation series itself. 
This problem has widely been discussed in the literature and many possible solutions have 
been proposed. Suggestions for the improvement of the analysis include the resummation 
of 7r 2 corrections using analytical continuation || (see also Q), the analysis of scheme 
dependence ||, and the inclusion of possible infrared renormalon contributions |J. 

In this letter we present an analysis along the lines of |3], |j] which consists of a partial 
resummation of the perturbation series for the semileptonic r width. We take advantage 
of the new result on the fourth order coefficient C3 of the beta function JF) in the series 
expansion 

(3(a) = -a 2 (l + ca + c 2 a 2 + c 3 a 3 + . . .) (1) 

where a = /5oa s /47r = 9a s /4ir for N c = 3 and Uf — 3. The specific choice of normalization 
of the coupling constant a used in Eq. ([I]) is convenient in higher order applications of 
the renormalization group since it removes all factors (3q from the relevant perturbative 
formulas. As an additional bonus the magnitude of all numbers become of order unity 
in the perturbative expansion. The connection with the traditional normalization is q = 

/V(/3o) m - 

The expression for the semileptonic r decay ratio has the form |§ 

Rt = r j r -* * + ha _ drons > = r- W r { -L) ^ (2) 



r(r — > v T + [i "*t 
with the weight factor 

W T (x) = 2(1 -x) 2 (l + 2x). (3) 

The spectral density p(s) = Jmil(s)/7r is defined in terms of the correlator H(Q 2 ) for 
charged semileptonic currents (see e.g. 0, 

We define a reduced r decay ratio r T by writing the semileptonic r decay ratio in a 
factorized form as 

R T = R° T (l + ^r T ). (4) 

R® is the lowest order partonic value of the ratio R T which is given by the number of 
colours N c to lowest order in the electroweak interaction. Using the experimental value 
Rt xp = 3.649 ± 0.014 |H| and = 3 one calculates 



T 



r^ xp = 0.487 ±0.011. (5) 
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It is well known that the reduced ratio r T in Eq. (^) has the power series expansion 

r T = a(m 2 ) + (6) 

We stress that corrections due to power suppressed terms, electroweak interactions 
and other possible corrections (mass terms, for instance) can easily be taken into account 
by a readjustment of R° T leaving the analysis of the hadronic part r T of the perturbation 
series untouched. A readjustment of R° T naturally implies a readjustment of the numerical 
value of r^ xp which is the only input parameter for our numerical estimates in massless 
QCD with n f = 3. 

We briefly outline the basic notions of the resummation analysis of the perturbative 
series of r T . The analysis is done within the MS scheme. Needless to say the resummation 
analysis can of course be performed in any other given scheme. In the resummation 
analysis one resums all those terms in the perturbation series that are accessible through 
the perturbative solution of the renormalization group equation || || . 

To start with, we define a reduced Adler's function d{Q 2 ) in analogy to Eq. (|J) by 
writing 

z,(<?2) = i( 1 + 5 <i(<?2) ) (7) 

for Adler's function D(Q 2 ). In the Euclidean domain this reduced function has the ex- 
pansion 

d(Q 2 ) = a(Q 2 ) + M 2 (Q 2 ) + ha 3 (Q 2 ) + k 3 a A (Q 2 ) + ... (8) 

The reduced rate function r T can then be expressed with the help of the /3-function 
coefficients c, q in Eq. (]!]) and the coefficients ki in Eq. (|j). One has |2j 
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where a T = a(m 2 ). In the expansion Eq. (|S]) we have explicitly listed terms up to and 
including a\. Up to 0{a%) all coefficients are known. For completeness we have also listed 
the a\ contribution with its unknown fc 3 -piece. Numerically one has 



a T + 2.312a 2 + 5.208a 3 T + (6.848 + k 3 )a* 



(10) 



where we have omitted the ellipses still present in Eq. (Q). For the numerical evaluation 
in Eq. flT0| ) we have used (see e.g. ||) 
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81 4374 



0.883 



'111 



3 



as well as 

h = ™- 4C(3) - 0.729 and k 2 = ^ - ^C(3) + ^C(5) - 1.258 (12) 

The five-loop coefficient k% is not yet known. Later on we present a simple estimate of k% 
which is used to obtain a feeling for the sensitivity of our results on the value of k^. It is 
noteworthy that the recently calculated four-loop coefficient |7] of the /^-function, 

140599 3560 > , , , . 

c 3 = + C 3 « 1.843 13 

118098 6561 SV ; V ; 

does not show up yet in this fixed order of perturbation theory (up to 0(af)). The 
coefficient C3 does, however, enter in our resummation analysis. Note that the series in 
the square brackets in each line of Eq. (|9|) depend only on the running of the coupling 
constant determined by the /3- function and the weight function W T in the integral (|2]). 
Therefore the corresponding "horizontal" series are accessible to resummation. This has 
been done before up to the third order in the /3-function ||, |4j]. In the present paper we 
go beyond the analysis in ||, |] by including the recently calculated four-loop coefficient 
of the /3-function [[7]. 

The results of resumming the series in Eq. ((^) can in general be presented in the form 

rv = a T M (a T ; (3) + k x a 2 T M x {a T ; (3) + k 2 a*M 2 (a T ] (3) + k z a A T M^a T ; /?) + ... (14) 

where the expansion coefficients Mi(a T ; (3) have no implicit dependence on the coefficients 
ki and reflect the effect of the running of the coupling constant and the dependence on 
the specific weight function W T (x) in Eq. (12). Explicit fixed order 0(a*) expansions of 
these coefficients are given in Eq. (||). The coefficient functions M$(a r ;/3) are completely 
determined by the /3-function and can be calculated explicitly using the renormalization 
group equation for the running coupling constant a with the starting value a T . Thus the 
coefficient functions can be obtained to all orders in the coupling constant a T in resummed 
form using the appropriate approximation for the /3-function. We emphasize that after 
resummation the formal ordering in terms of powers of the coupling constant a T in Eq. (|S|) 
is lost (compare [|K|). 

Let us return to the integral representation Eq. (fj) which we shall rewrite in reduced 
form according to Eq. (|j). Accordingly we define a reduced part of the correlator by 
p(Q 2 ) and integrate p(Q 2 ) along the contour \s\ = m 2 . in the complex Q 2 = —s plane. We 
obtain 

WT ( 1 \rJ S )± = ir W JJL)^ ( 15 ) 



The function p{— Q 2 ) is calculated by integrating the differential equation 



d(Q 2 ) = -Q 2 ^^ 1 - (16) 



Some words of caution are in order concerning the use of Eq. ([15]). In the second part 
of Eq. flT5|) we have assumed that one can use the perturbative approximation for p{Q 2 ) 
along the contour \s\ = m 2 . Due to the existence of Landau singularities (corresponding 
to poles even in lowest order of perturbation theory) the perturbative expression for p(s) 
does not have the correct analyticity properties as prescribed by the Kallen-Lehmann 
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representation of the full correlator function. The result of doing the contour integral 
therefore also depends on whether the contour encircles the Landau singularity or not. 
An alternative definition of the integral on the left hand side of Eq. (|i~5"D would be to use 
the discontinuity of the renormalization group improved correlator function p(Q 2 ) across 
the positive semi-axis. These two definitions of the integral of the renormalization group 
improved spectral density differ by an amount (A 2 /m 2 ) n , where A is the usual scale of 
QCD and n is determined by the power of x in the weight function Eq. (||). In perturbation 
theory one has a pure polynomial expansion in terms of the coupling constant. As a result 
of analytically continuing the polynomial terms in the complex plane one obtains factors 
of (\nz) n . They all have the correct analyticity properties (cut along the positive semi- 
axis) and the above ambiguity does not occur in any finite order of perturbation theory. 
Redefinitions of this kind (or direct treatment of running on the positive semi-axis) were 
considered in the literature |TI|, [L2|, [L3|. We do not further dwell upon this point but 
consider the right hand side of Eq. (|15"D as a definition of the left hand side keeping in 
mind the above possibility of an alternative definition. 

In order to proceed with the resummation of the perturbation series (for early refer- 



ences see e.g. |]I1|) we decompose the coefficient functions Mj(a r ; (3) into moment functions 
according to 

n + 1 f , ,„ 2 



M hn (a T ; (3) = — — i {-x) n a\m 2 x)dx. (17) 
2ixi J\x\=i 

Considering the weight function W T in Eq. (|3]) one has 

M,(a r ; (3) = 2M ifi (a T ; (3) - 2M !i2 (a T ; (3) + M ij3 (a T ; (3). (18) 

These moments constitute the building blocks of our analysis. 

As an instructive example we consider the resummation of the moment function 
Mi <0 (a T ; (3) to leading order in the /3-function, i.e. f3(a) = —a 2 . To this order of the 
/3-function, the moment M 1)0 (a T ;/3) corresponds to the a 2 (Q 2 )-term in d(Q 2 ) which in 
turn is contributed to by the a(Q 2 )-term in p(Q 2 ) such that 



<M li0 (a T ; (3) = - — i a(-m 2 x)dx = — / \ . (19) 

2m J\x\=i Z7T J-tt 1 + ia T (p 

One finally has 

M li0 (a T ;/3) = J- f (20) 
27ca T J-tt 1 + ia T <p 

where a T = a(m 2 ) = (ln(m 2 /A 2 )) -1 is the solution of the renormalization group equation 
to lowest order. The generalization of this technique to higher orders in the /3-function 
and the reduced Adler's function d{Q 2 ) in Eq. ([8]) is quite straightforward 0, fl| and 
will not be described in much detail. We have extensively used MATHEMATICA [[L5j 
for the fourth order analysis and found it to be a very appropriate tool for doing these 
calculations. 

In this way we can determine values of the strong coupling constant using the re- 
summed form of Eq. (^) as written down in Eq. ( |H| ) . In order to exhibit the dependence 
of the resulting strong coupling constant on the order of the /3-function accuracy we list 
the resummed result for increasing orders of /3-function accuracy. We obtain 

4 1} = 0.2666 ±0.0041, a {2) = 0.2668 ± 0.0045, 

4 3) = 0.2687 ±0.0048, of = 0.2704 ± 0.0050 (21) 
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or, in terms of a^jg, 

a® = 0.3723 ± 0.0058, af ] = 0.3726 ± 0.0063, 

4 3) = 0.3751 ±0.0067, a[ 4) = 0.3775 ± 0.0071 (22) 

where refers to the use of the first order /3-function (3(a) = —a 2 , etc. For the sake 
of comparison we also give the result of calculating a T to fixed order in perturbation 
theory using Eq. ([10|) to third order in a T . One obtains a^y = 0.2535 ± 0.0035 (i.e. 
a® = 0.3540 ± 0.0048). When we quote errors on these numbers the error is entirely 
determined by the experimental error on r£ xp as given in Eq. 

We postpone the discussion of the numerical results in Eqs. ([H]) and (p2]) until after 
analyzing the convergence properties of the resummed functions. The first attempt to 
determine the convergence radius of the resummed functions has been done in ||. The 
authors of [Q obtained approximate values for the convergence radius by numerical inves- 
tigation. Here we present a new method to determine the convergence radius by analyzing 
the singularity structure in the complex a T -plane. 

In the simple lowest order example given above by Eq. (|T9|) , the radius of convergence 
is determined by the solution of the equation 1 — ma = 0, leading to the region of 
convergence \a\ < 1/tt 0. In higher orders of /5-function accuracy the investigation of 
the convergence properties of the resummed functions M i n (a T ,/3) in Eq. (O) is not as 
simple. It is clear that the numerical investigation employed by the authors of |4j] is no 
longer efficient as one gets closer and closer to the true radius of convergence. In fact one 
can do better than that. The key observation is that we can search for convergence radii 
as singularities of the resummed function M^ n {a T , (3) continuing a T to the whole complex 
plane. The evolution of a along the contour Q 2 = m^e 1 ^ (0 G [— 7r, 7r]) is governed by the 
renormalization group equation 

da 

^ttt = Pip) = ~a 2 (l + ca + c 2 a 2 + c 3 a 3 + . . . ) (23) 

U(p 

with any given complex a T as starting value. When a T is real, it is clear that there is no 
singularity of the solution of the renormalization group equation ( ]2"3"| ) and thus no singu- 
larity of the resummed functions M^ n {a T1 0). In fact, the search for singularity solutions 
of Eq. Q2"5p has to be extended to the whole complex a T -plane. The closest singularity in 
the complex a r -plane then determines the convergence radius of the resummed functions 
M iin (a T ,/3). 

Technically one searches for the singular solutions of the renormalization group equa- 
tion (^) by analyzing the solution to ( p3|) in implicit form, 

i(j)(a;a T ) = f 

Ja T p[a') 

The integral is a line integral in the complex a-plane connecting the complex starting 
value a T with its evolved value a. Eq. ( |2"4| ) determines the value of the angle in the 
evolution from a T to a with 0(a = a T ; a T ) = 0. Singularities in the evolution process are 
encountered at angles (p s (a T ) = 0(a = oo; a T ). After inversion the singular starting values 
a T ((f) s ) appear as lines in the complex a T -plane which are parametrized by the angles <p s at 
which the singularity a = oo occurs. The search for the singularity lines in the complex 
a T -plane can be conducted with the help of built-in facilities of MATHEMATICA. For 
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Table V. Higher order perturbative predictions for the reduced part of the semileptonic r 
decay width using a^> = 0.2535 



example, in Fig. 1 we have drawn the singularity lines a^((j) s ) in the a r -plane for four- 
loop /^-function accuracy. The singularity closest to the origin determines the radius of 
convergence of the resummed series Mj in (a T ) to four-loop /^-function accuracy. 

Let us now list the values of the radii of convergence for increasing orders of /3-function 
accuracy including the one-loop result. One has 



a« = 1/tt = 


0.318 




= 0.444), 


4 2) = 0.744/tt = 


0.237 


(4 2) 


= 0.331), 


4 3) = 0.697/tt = 


0.222 


(4 3) 


= 0.310), 


4 4 ) = 0.674/tt = 


0.214 


(4 4) 


= 0.299), 



where the quoted numbers for d@\ and a)y have been calculated up to three-digit 
accuracy. We mention that the algorithmic procedure allows one to determine the radii of 
convergence to any desired degree of accuracy. We have added the corresponding numbers 
for the conventional coupling constant a T = a s {m^) = 4ira T /9 in brackets. The authors 
of [|J have determined an upper bound for the three-loop convergence radius a£ ■ They 
quote 4- < 0-25 which is 12% above the true convergence radius a® listed in Eq. (p5|) . 

Looking at the numbers in Eq. (p5|) one sees that the convergence radii become smaller 
as the order of the /^-function accuracy increases. With the newly calculated four-loop (3- 
coefficient C3 the convergence radius is reduced by 3.6% relative to the three-loop result. 
The 3.6% reduction in the convergence radius must be contrasted with the mere 0.6% 
increase in the value of a T when going from three-loop to four-loop order in the /3-function 
(see Eq. (^)). It is interesting to speculate that the shrinking of the convergence radius 
continues as one goes to ever higher orders of /3-function accuracy including the possibility 
that the convergence radius shrinks to zero when the order of the perturbative /^-function 
expansion goes to infinity. 

As mentioned after Eq. (|22|) , the fixed third order value of a T is given by af> = 
0.2535 which clearly is outside the convergence region. This means that the perturbative 
approximation for moments diverges at the scale determined by the experimental data 
for the semileptonic r decay width. The resummed results for a T in Eq. (ET]) lie outside 
the respective convergence regions. The message is clear. The resummed values are not 
accessible by using higher and higher order approximations of perturbation theory. 

Let us illustrate the lack of convergence for the fixed order perturbation series of r T 
using this value a® = 0.2535. In Table [l] we list the resulting values of r T for increasing 
orders in a T for this case. The n = 3 value for r T is fixed by the experimental input 
value r^ xp = 0.487 by construction. Nothing dramatic happens to the series up to n = 8. 
Starting with n = 9 the series can be seen to show a divergent behaviour. 
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Table 2: Resummed moments for f3(a) = — a 2 (l + ca + c 2 a 2 ) 
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Table 3: Resummed moments for (3{a) = — a 2 (l + ca + c 2 a 2 + c 3 a 3 ) 



Let us now return to discuss the resummed coupling values in Eq. ([H]). Increasing 
the /9-function accuracy does not appear to change the numerical results much. In order 
to exhibit some of the structure of our resummation analysis we present results on the 
moment functions Mj >n (a; f3) and their relevant linear combination Eq. (|T8|) in Tables 
and In Table |2| we have used three-loop /3-function accuracy relevant for the determi- 
nation of the resummed value of in Eq. (^1|) while Table |3| contains our result on the 
resummed value of up to four-loop /3-function accuracy. For the sake of comparison we 
list in Table f| values for the fixed third order expansion for the corresponding moments. 
Comparing the numbers in Tables |2| and |3| one sees a weak dependence on which order of 
the /3-function is used. The dependence is sowewhat larger for the coefficients i = 1 and 
i = 2, but the contributions of these coefficients to the total result are suppressed. The 
bulk contribution comes from the leading order term (i — 0). Comparing Tables |2] and ^] 
on the one hand and Table [| on the other hand one sees that the effect of resummation 
is largest for the higher order coefficients i = 1,2. However, the effect of the higher order 
coefficients on the total result is small. The overall effect of resummation at this level is 
not significant. Resummation does, however, lead to a slight increase of the value of the 
strong coupling constant extracted from r decay data though, namely from a s = 0.3540 
to a s = 0.3775. Nevertheless, such an increase must be taken seriously considering the 
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Table 4: Moments for the fixed third order expansion 
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level of accuracy that the ct s determination has reached in other experiments. We shall 
return to this point later on. 

The value of k 3 is still unknown. One can obtain a first estimate for the value of k 3 by 
means of a Pade approximation (e.g. [JT6|). Using this estimate one can get a hint on the 
sensitivity of our results on higher order corrections in the MS scheme. The naive (1, 1)- 
Pade approximation for the function d{Q 2 ) results in the value k 3 = k\jk\ = 2.17 [|, |T7| . 
Performing the resummation analysis with this value of k 3 one finds a value a^(k 3 = 
kl/h) = 0.2650 (a^(k 3 = k 2 /h) = 0.3700). 

We want to emphasize that the renormalization group equation and therefore our 
resummation procedure is invariant under a one-parameter subgroup parametrized by 7 
which describes the rescaling of the renormalization scale m T . The coupling parameters 
in the two schemes are related by 

a\m 2 T ) = a{e<m 2 T ) = a - 7a 2 + ( 7 2 - c 7 )a 3 - ( 7 3 - ^c 7 2 + c 27 )a 4 + . . . . (26) 

Rewriting the reduced d(Q 2 ) function in terms of a' leads to a rescaling of the coefficients 
ki but leaves the /3-function invariant. A suitable choice of the parameter 7, which cor- 
responds to a different choice of renormalization scheme, may lead to a value a' for the 
coupling inside the convergence region. At the same time one may have deteriorated the 
series in ki. This possibility underlines the fact that the coupling constant in a particular 
scheme has no physical meaning. 

To get rid of this ambiguity, we choose a physical observable, namely the reduced 
Adler's function itself as an effective coupling constant [IB, |T9"|| , 



d(Q 2 )=a k (Q 2 ). (27) 
For this effective coupling constant the /3-function is given by 

Pk(a k ) = -a 2 k (l + pidfc + p 2 a\ + p 3 a 3 k ) (28) 

with 

Pi = c, p 2 = c 2 - cki -k\ + k 2 , 

p 3 = c 3 - 2c 2 h + ck\ + 4A; 3 - 6hk 2 + 2k 3 . (29) 

Because the /3-function is expressed in terms of a physical observable, the coefficients pi, 
p 2 and p 3 are renormalization group invariants [^] (see e.g. |2l|). Now there is no gain in 
knowing c 3 in the MS scheme because the four-loop coefficient c 3 of the /3-function appears 
together with the unknown coefficient k 3 in the invariant p 3 . This observation reflects the 
fact that one needs full information on the q and ki for a given physical observable in any 
given renormalization scheme. Numerically we obtain p\ m 0.790, p 2 ~ 1.035 and 

37096148 4820288 ^ N 12352 >/ x9 > , N o 59800 > . , 

p 3 = + C(3) + C 3 2 - 256C3 3 C(5) 

H 59049 6561 v ; 81 w w 243 v ; 

U>00 -C(3)C(5) + 2k 3 = -2.97953 + 2k 3 w 1.360 (30) 
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using the Pade- approximation k 3 = k\jk\ = 2.17. The region of convergence in this 
scheme can be calculated to be |a^| < 0.214 and again does not include the value = 
0.3630. 
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From the explicit results in Eq. ( ]50j ) we also see that there is an interesting possibility 
of having the infrared fixed point in the Euclidean region if the coefficient k% happens to 
be small enough. In naive perturbation theory the sign of k$ is not very important (as well 
as in the resummation approach within the MS scheme) while in the effective coupling 
scheme things can change qualitatively for negative p%. Because physics is determined by 
invariant quantities this would mean that the scheme dependence stops to be soft in this 
order of perturbation theory. 

Let us summarize our main numerical results. We have analyzed semileptonic r decay 
data using a partially resummed perturbation series and have obtained a value of the 
strong coupling constant which is approximately 7% higher than the value obtained from 
the fixed third order analysis. In the resummed analysis we find a-g^ = 0.3775 compared 
to a^jg = 0.3540 in the fixed order expansion analysis. When we rescale the results up 
to the Z mass, the resummed value of from r decay tends to be in better accord 
with the determination of a s from Z decay data |J. This is good. On the other hand a 
higher value of a s is not welcome when one does low energy sum rule phenomenology 



Even the present fixed order value a s = 0.3540 is rather large when one wants to justify 
the perturbative computation of coefficient functions in the framework of the operator 
product expansion. 

Judging from the results presented in this paper there does not seem to be much 
room for improvements on the accuracy of our results. Adding the new four-loop term to 
the /^-function has not changed the results of the resummation analysis by a significant 
amount (a^ compared to in Eq. fl2T|)). Apart from the fact that the five-loop result is 
difficult to come by we do not anticipate that adding the five-loop term to the /3-function 
will affect our resummation results much. We have attempted to estimate the influence of 
the five-loop k 3 term. The analysis has shown that higher order /^-contributions tend to 
be numerically suppressed. A pattern of such a suppression can be seen in Eq. (|IT]) with 
^3 = k\jk\ = 2.17. Although this is not a renormalization group invariant statement, 
there is the temptation to skip the higher order ki contributions altogether and use the 
resummed remainder as an improvement of perturbation theory. 

It is very likely that our approach makes more sense than to push fixed order pertur- 
bation theory to ever higher orders. For example, the fixed fifth order formulas that have 
been written down in 17] within the latter approach may not even lead to an improved 
analysis. As our numerical results have shown, the resummed values of a s from r decay 
lie outside the convergence radii and can therefore not be obtained from a power series 
expansion. Regular perturbation series do not converge to their resummed counterparts. 
The experimental value of r T appears to be too large for a fixed order perturbation analysis 
to apply. 

We end our discussion with some critical remarks. The simple analyticity properties 
of the two-point correlator (Kallen-Lehmann representation) describing semileptonic r 
decays are central to the present resummation analysis. In other more complicated cases 
a corresponding resummation analysis may not even exist. We do not even attempt to 
answer the question whether it makes sense to compare two values of a s derived from 
different experiments, where the one value is extracted using a resummation analysis and 
the other value is obtained using fixed order perturbation theory. 

Last but not least, we reiterate that the entire scheme presented here is based on a par- 
ticular prescription of analytic continuation in the complex Q 2 -plane which corresponds 
to a particular prescription for treating the Landau pole inside the integration contour. If 
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one chooses another prescription, the derived values for will change. This uncertainty 
is inherent to the perturbative solution of the renormalization group equation since the 
perturbative solution has wrong analyticity properties in the complex Q 2 -plane. 
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Figure Caption 



Fig. 1: Singular starting values a^((f) s ) ((f) s E [— it, 7r]) as parametric curves in the 
complex plane. The radius of the dashed circle determines the convergence 
radius of the resummed series. The resummed value = 0.2704 extracted 
from the r decay rate lies outside the region of convergence. 
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